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Abstract 

We consider the decay of an intermediate mass Higgs boson into hadrons up to order 
O(a^). The results from the diagrams containing only light degrees of freedom were 
recently computed in a previous work by K.G. Ch. (Phys. Lett. B 390 (1997) 309). 
In this letter the remaining contributions involving the top quark are evaluated 
analytically using an effective field theory approach. Coupled with the previous 
result the present calculation determines the complete next-next-to-leading order 
correction to the hadronic decay rate of an intermediate mass Higgs boson. 


The Higgs boson is the only still missing particle in the standard model of particle 
physics. Up to now only mass limits exist from the failure of experiments at the CERN 
Large Electron-Positron Collider (LEP1) to observe the process e + e~ —> ZH. Currently 
the mass range Mh < 65.6 GeV is ruled out at the 95% confidence level [|T[ . 

From the theoretical side a lot of effort has been invested to get more insight into the 
properties of the Higgs boson (for a review see j2|). Of particular interest is thereby the 
intermediate mass range where mt, -C MR ^ 2 M\y. Then the dominant decay mode is 
the one into bottom quarks. Concerning QCD corrections the full mass dependence at 
0(a s ) was evaluated in 0. However, especially for the fermionic decay of an intermediate 
mass Higgs boson it turns out that an expansion in the quark masses provides reasonable 
approximations. At O(o%) besides the massless limit (i.e. keeping only the overall factor 
m b) 0 a ^ so sub leading mass corrections in the rril/Mjj expansion are known [0-0- 

Recently the contribution to the scalar polarization function of the Higgs boson con¬ 
taining only light quarks was considered at four loops 0. Its imaginary part was evaluated 
analytically in the massless limit leading to corrections of O(a^) to the Higgs decay rate. 
The corresponding corrections prove to be numerically more important than the power 
suppressed contribution of O{ot 2 s ml / Mff). 
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In Ref. [6] additional quasi-massless (and numerically important) contributions of order 
a 2 have been identified and elaborated. They come from so-called singlet diagrams with 
a non-decoupling top quark loop insideQ The results of Ref. || have been confirmed 
and extended by the computation of power-suppressed terms of order a 2 s (Mjj / M 2 ) n (n = 
1, 2,...) [9;]. For completeness we should mention that the exact result for the imaginary 
part of the double-bubble diagram with massless external quark and internal top quark 
can be found in |H|]. Radiative corrections enhanced by a factor M 2 , usually expressed 
through X t = GpM 2 /8 tt 2 \/2, are also available up to the three-loop order [11], [12| and 
will be compared at the end with the new terms of O(a^). 

In this letter we calculate the top-induced corrections at 0(a^) and hence complete the 
analysis of the hadronic Higgs decay at NNNLO, as far as one neglects power suppressed 
corrections. Although the phenomenological interesting decay is the one into bottom 
quarks the following discussion will be kept more general and a generic light quark q will 
be considered. However, for the numerical discussion we will come back to the case of 
bottom quarks. 

We start with the bare Yukawa Lagrangian, 


jC y = — ( y^m°q°q 0 + 


( 1 ) 


where v is the Higgs vacuum-expectation value and the superscript 0 labels bare quantities. 
Assuming that the Higgs boson mass, Mh is less than the top quark mass m t , Cy can 
be replaced by an effective Lagrangian produced by integrating out the top quark field. 
According to Refs. 


13, 


the resulting Lagrangian reads 


/^eff _ 

— 


H°_ 

n° 


Cl [Oil + £ c 2q [O', 


( 2 ) 


where [0\] and [0' 2q ] are the renormalized counterparts of the bare operators 


o; = (g®„) , o', = m®yy>', 


with v being the (bare) field strength tensor of the gluon. The primes mark the 
quantities defined in the effective nq — 5 QCD including only light (in comparison to the 
top quark) u, d, s, c and b quarks. All the dependence on the top quark gets localized in 
the coefficient functions C\ and C^q■ If the latter are known the computation of the total 
decay rate Th = T(H —> hadrons) is reduced to the evaluation of the functions A jk(M^) 
(jk = 11,12,22) related via the optical theorem to the absorptive parts of the scalar 


correlators 

n Al 2 ) = i J dxe iqx (0\ T\ 0^)0^ 0) ] |0> 


(3) 


It should be stressed that by the very meaning of the effective Lagrangian (j2j) the cor¬ 
relators of Eq. ([|) may be computed within the effective massless QCD, which leads to 

1 Similar effects for the decay of the Z-boson to quarks were first discovered in J^]. 
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Figure 1: Typical Feynman diagrams contributing to A 22 . 
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The solid circles represent 




Figure 2: Typical Feynman diagrams contributing to The solid circles represent 

the operator [0(]. 


a drastic simplification of the calculation. In addition, the coefficient functions C\ and 
C 2 q essentially depend on only one kinematical variable, the top quark mass, which also 
simplifies their calculation a lot. 

Typical diagrams contributing to A 22 , An and Af 2 are shown in Figs. [l|-[| Note that 
every massless diagram contributing to A 22 to order af does only have cuts containing 
at least one quark-antiquark pair. This allows one to unambiguously associate A 22 to 
T(H —> qq) — the partial decay rate of the Higgs boson to hadrons containing at least one 
qq pair. The diagrams contributing to A n (see Fig. §) describe the production of gluons. 
The leading order diagram has clearly only contributions to the gg final state. Starting 
from NLO, however, there are contributions from diagrams leading both to gg, ggg and 
gqq final states. The interpretation of the Af 2 (see Fig. |3|) is even more complicated. 
Here different final states begin to appear already in the leading non-vanishing order |6|. 
In this paper we will adopt a pragmatical point of view and evaluate the total hadronic 


decay rate without any attempt to differentiate between final states [15 


The functions An, C\ and A 22 have been computed analytically in |l3], |l(| [14| and 
0. Thus, to complete the evaluation of Th in the O(o%) approximation one needs to 
compute C 2( j as well as A^ 2 . This calculation and its results are described below. 

For convenience we should recall the relevant formulae which connect the bare and 
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Figure 3: . Two- and some of the three-loop diagrams contributing to Af 2 - The solid 
circles represent the operators [0(] and 0 2q , respectively. 


renormalized quantities appearing in Eq. (Q). For the operators we have 


m = 


1 + 2 


' a s d 

da. 


In Z n 


O' = O' 


'2q 


2 qi 


' a, d 


0',- 4 (^ 1 nZ m )E+. 


(4) 


and for the coefficient function the relations look like: 


Q = 


C 2q = 


■Cl 


1 + 2 (a s d/da s ) In Z, 
A(a s d/da s )\wZ m 0 0 

l + 2(a s d/da s )\nZ g 1 2q ’ 


(5) 


where the renormalization constant Z g and Z m are defined in the effective theory. In 
the order we are interested in these constants are needed up the one- and two-loop level, 
respectively, because is already proportional to a s . 

For our purpose we need C\ up to order a 2 which may be found in JO], |H| [Hj . C 2q is 
also known up to O(a^) [|1^] but for the present analysis needed up to 0(a,). To this aim 
we want to use the low energy theorem (LET) which allows to attach the Higgs boson 
via differentiation w.r.t. the masses involved in the process. The formula given in |12| for 
the computation of C 2q simplifies in our case to 


q 


= 1 


g(+(°) + £ 


( 0 )). 


( 6 ) 


It is understood that after the derivative w.r.t. the bare top mass is done the renormaliza¬ 
tion of the parameters m t and a s is performed. The superscript t indicates that only the 
diagrams where at least one top quark is present have to be taken into account. The first 
non-vanishing contribution arises at two-loop level involving one diagram. At three loops 
altogether 25 diagrams contribute. Some typical examples are depicted in Fig. £f|. Note 
that only the pole parts have to be computed because the integrals are logarithmically 
divergent and consequently proportional to (/i 2 /mj) £ . The differentiation w.r.t. the top 
mass then leads to an additional factor e. The diagrams are generated with the program 
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Figure 4: Two- and some of the three-loop diagrams contributing to Eg and E y. The 
thick (thin) lines represents the top quark (light quarks). 


QGRAF JT7| and fed into the package MATAD written in FORM for the purpose to 
solve three-loop tadpole integrals. The calculation was performed with arbitrary QCD 
gauge parameter £. Eg and E^ separately still depend on £, however, the final result for 
C* 2 q is independent of the gauge parameter which serves as a welcome check. Expressing 
the result in terms of the MS top mass, rrit, we get 


C 2q — 1 + 


«i 6 ) (h)' 


7T 


---ln^- 

18 3 m{ 


+ 




7T 


311 

Ll296 + 3^ ^ 


175, u 

-In — 

108 


2 29 2 /j 2 

2 _ FT ln —2 + n i 
mf 36 


53 1 9 u 

- 1 -In 2 — 

216 18 


m? 


(7) 


with £(3) ~ 1.20206. rq is the number of light quarks. Using the relation between the 
MS and the on-shell mass, M t , at one-loop level, mt/M t = 1 — (4/3 + In /i 2 /M 2 )a s / 7 r, one 
gets: 


C? q S = 1 + 


v 71 . 

247 


—-- In -^-2 

18 3 M 2 


+ 


, /i 2 29, , n 2 

In —^-In" —k + ni 

108 M 2 36 M 2 




7r 


111 + -C(3) 
1296 3 v 2 


216 18 M 2 


( 8 ) 


For completeness we also list the result for C\ 


71 


c = 1 af M 

1 12 IT 


1 + ^(H_l ln 4' 

7r l 4 6 m 2 


(9) 


On the basis of the effective Lagrangian of Eq. (D) it is possible to write down the 
decay rate of the Higgs boson into hadrons in the following form: 


r (H - hadrons) = (1 + S u ) 2 I A qq [(1 + A^ 2 ) (C 2q ) 2 + Af 2 C' 1 C' 2g + A 


hdo 


+ Agg AH (Cl)" 


( 10 ) 
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with A q q = 3 GfM/Att\/2 and A gg = AGfMH'k\/2. The universal corrections S u arise 
from the renormalization of the factor H°/v°. The terms up to 0(a 2 X t ) can be found in 
11], A hdo summarizes the corrections coming from higher dimensional operators. They 
are at least suppressed by a factor a 2 * Ml/M 2 . The factors A tJ contain the electroweak 
and QCD corrections from the light degrees of freedom only. As we are interested up 
to an accuracy of 0(a^) An and A 22 are known. However, A\ 2 is only available at the 
two-loop level where one diagram has to be evaluated (see Fig. §). We have to extend 
the analysis to three loops where altogether 28 massless three-loop diagrams have to be 
computed. This was done with the help of the program package MINCER |0|. In Fig. [| 
some graphs are pictured. The single diagrams again depend on £ which cancels in the 
proper sum. The result reads: 


Af 2 = 


« 1 5 ) (r) 


7 r 


“_ 8 ln-£. 
3 Ml 


H 


+ 




7T 


l0 ° 73 + 76C(2) + 156C(3) - In -Ar - 38 In 


2 I 1 


18 


M 2 H 


M 2 h 


(283 8 . . 16 . . 112 fj, 2 

+ m -C(2)-C(3) H-In -A— 

9 3 SV ; 3 SV ; 9 Ml 


H 


+ -ln 2 ^4 

3 M 2 H/ 


( 11 ) 


with £(2) = 7r 2 /6. The 0(a s ) terms can also be found in [[], j9L ]T2[ , the 0(a 2 ) terms are 


new. 


Let us now compare the new results of O(a^) with the previous ones and also with 
other known correction terms which might be of the same order of magnitude. Here, 
we have in mind terms of order Oi 2 rn b / Mfj , aa s , oi 2 X t and a 2 s M'l/M 2 . Thereby it is 
convenient to write Eq. (|T0|) in the form 2 : 


F (H — > hadrons) — A bb ^1 + A^ + A t ^ + (1 + S u ) A gg An (Ci) , 


( 12 ) 


where A b t contains only corrections form light degrees of freedom and all top-induced 
terms from the first line of Eq. ([II]) are contained in A t . Furthermore we express A t also 
in terms of q1 5 )(/t). Choosing fi 2 = Ml and ni = 5 we find 


A" - _ 6 MT 

' “ ® Ml 


+ 0.472 Aisl + 0.651 AAl a g> + a g) | 5.667 - 40.000 lm '’ 


( 5 )\ 2 ' 


7T 


7T 


( C 


{ S) 


29.147-87.725 




( 5 )\ 2 ' 


Ml 


41.758 ( 


A t = I 3.111 — 0.667 L t + 


(m 


( 5 ) 4 2 


4 ] t 

4 


Ml 


(13) 


Ml 


— 10 + 4 L t -\— In 



2 From now we will consider all quarks with mass lighter than m b as massless. This means that the 

sum in Eq. ( 0 ) reduces to q = b. Furthermore the numerical discussion comparing the new corrections 

with previously known terms is restricted to the part proportional to to A bb . 
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( 14 ) 


3 2 M 2 

+ (ag } ) (50.474 - 8.167 L t - 1.278 L 2 ) + (ag } ) -^f (0.241 - 0.070 L t ) 
+ W (l-4.913ag } + (ag } ) 2 (-72.117- 20.945 , 


The term in Eq. (|b|) proportional to 

( mJ 5) ' 2/ »' 2 


with = cx^(Mh)/tt and L t = In Mfj/Mf. 

Mfj/M 2 is the leading contribution from A hdo . The (m y b 0> ) 2 /M'jj corrections in A t arise 
from the singlet diagram with one top and one bottom quark triangle and can be found 
in [|J. At this point we should mention that still contributions with pure gluonic final 
states are contained in Eq. At O(al), however, these corrections are not yet known 
and can thus not be subtracted. 

In the approximation considered in this paper we have —2 <, L t < 0. This means that 
the logarithm needs not necessarily to be resummed as in addition the coefficients in front 
of L t are much smaller than the constant term. 

One observes that the new top-induced corrections of 0(a d ) are numerically of the 
same size as the previous ones arising from “pure” QCD. Furthermore one should mention 
that the coefficient of the M t -suppressed terms are tiny and, as a s /X t ~ 30, also the a 2 X t 
enhanced terms are less important than the cubic QCD corrections. For comparison in 
Eq. flKf ) also the two-loop corrections of order aa s are listed. In principle also higher 
order mass corrections are available However, in the case of bottom quarks it turns 
out that they are tiny. 

To summarize, in this letter the top-induced corrections of order for the hadronic 
Higgs boson decay were presented. An effective Lagrangian was constructed and both 
the coefficient functions and relevant correlators were evaluated up to three loops. The 
new results combined with the ones in @ lead to complete corrections of 0(a d ) to the 
hadronic Higgs decay. 
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